/*=====================================================================*
 *                   Copyright (C) 2011 Paul Mineiro                   *
 * All rights reserved.                                                *
 *                                                                     *
 * Redistribution and use in source and binary forms, with             *
 * or without modification, are permitted provided that the            *
 * following conditions are met:                                       *
 *                                                                     *
 *     * Redistributions of source code must retain the                *
 *     above copyright notice, this list of conditions and             *
 *     the following disclaimer.                                       *
 *                                                                     *
 *     * Redistributions in binary form must reproduce the             *
 *     above copyright notice, this list of conditions and             *
 *     the following disclaimer in the documentation and/or            *
 *     other materials provided with the distribution.                 *
 *                                                                     *
 *     * Neither the name of Paul Mineiro nor the names                *
 *     of other contributors may be used to endorse or promote         *
 *     products derived from this software without specific            *
 *     prior written permission.                                       *
 *                                                                     *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND              *
 * CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES,         *
 * INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES               *
 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE             *
 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER               *
 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,                 *
 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES            *
 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE           *
 * GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR                *
 * BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF          *
 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT           *
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY              *
 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE             *
 * POSSIBILITY OF SUCH DAMAGE.                                         *
 *                                                                     *
 * Contact: Paul Mineiro <paul@mineiro.com>                            *
 *=====================================================================*/

#include <pthread.h>

#include "qfunc.h"
#include "prating.c"

float
qfunc (const float*                     loga,
       unsigned int                     rank,
       float                            logbeta,
       float*                           dlogbeta,
       float*                           ddlogbeta,
       const float*                     c,
       unsigned int                     n_labels,
       const float*                     d,
       const float*                     pz,
       const Rating*                    ratings,
       unsigned int                     n_ratings,
       const Distribution*              priorlogbeta,
       float                            hq)
{
  float q = hq;

  *dlogbeta = 0;
  *ddlogbeta = 0;

  for (unsigned int i = 0; i < n_ratings; ++i)
    {
      for (unsigned int j = 0; j < n_labels; ++j)
        {
          float dlogbetatmp = 0;
          float ddlogbetatmp = 0;

#define logafunc(n) loga[ratings[i].worker * rank + n]
#define cfunc(n) (1.0 + c[n_labels * n + j])
#define dfunc(n, k) ((k == j) ? 0 : (1.0 + d[n_labels * n + k]))

          q += pz[j] * 
               logprating (ratings[i].rating,
                           logafunc,
                           rank,
                           logbeta,
                           dlogbetatmp,
                           ddlogbetatmp,
                           cfunc,
                           dfunc,
                           n_labels);

          *dlogbeta += pz[j] * dlogbetatmp;
          *ddlogbeta += pz[j] * ddlogbetatmp;
#undef dfunc
#undef cfunc
#undef logafunc

        }
    }

  q += distribution_logpdf (priorlogbeta, logbeta);
  *dlogbeta += distribution_dlogpdf (priorlogbeta, logbeta);
  *ddlogbeta += distribution_ddlogpdf (priorlogbeta, logbeta);

  return q;
}

void
dqfunc (float*                          dloga,
        float*                          dc,
        float*                          dd,
        const float*                    loga,
        unsigned int                    rank,
        float                           logbeta,
        const float*                    c,
        unsigned int                    n_labels,
        const float*                    d,
        const float*                    pz,
        const Rating*                   ratings,
        unsigned int                    n_ratings)
{
  for (unsigned int i = 0; i < n_ratings; ++i)
    {
      for (unsigned int j = 0; j < n_labels; ++j)
        {
          float tmpdd;

#define logafunc(n) loga[ratings[i].worker * rank + n]
#define dlogafunc(n) dloga[ratings[i].worker * rank + n]
#define cfunc(n) (1.0 + c[n_labels * n + j])
#define dcfunc(n) dc[n_labels * n + j]
#define dfunc(n, k) ((k == j) ? 0 : (1.0 + d[n_labels * n + k]))
#define ddfunc(n, k) *((k == j) ? &tmpdd : dd + n_labels * n + k)

          dlogpratingdall (pz[j],
                           dlogafunc,
                           dcfunc,
                           ddfunc,
                           ratings[i].rating,
                           logafunc,
                           rank,
                           logbeta,
                           cfunc,
                           dfunc,
                           n_labels);
        }
#undef ddfunc
#undef dfunc
#undef dcfunc
#undef cfunc
#undef dlogafunc
#undef logafunc
    }
}
